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ABSTRACT: The original Jeans dispersion relation and instability criterion are derived 
by a mathematically well-defined limiting procedure. The procedure highlights Jeans' 



t^ I physical reasoning and vindicates the (in)famous "Jeans swindle." A second, independent 

v" ■ procedure is stated which yields the same result. 
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1. INTRODUCTION 

In 1902, J.H. Jeans [1] derived his celebrated dispersion relation 

uj^ = \kfcl-4nGpo (1.1) 

governing the evolution of infinitesimal disturbances of a fictitious infinitely extended, 
homogeneous and isotropic, static fiuid of mass density po that is coupled to Newtonian 
gravity. According to (1.1), an initial disturbance whose wavelength 27r/|k| is much smaller 
than the Jeans length 

behaves essentially like a classical sound wave of sound speed Cg, with self-gravity con- 
tributing only slight corrections, but, as Jeans noted [1], perturbations whose wavelengths 
surpass the Jeans length can grow exponentially in time. The dispersion relation (1.1) has 
a counterpart in the kinetic theory of encounter-less stellar dynamics; see, for instance, 
the monograph by Fridman and Polyachenko [2]. A related yet conceptually somewhat 
different question is the evolution of infinitesimal disturbances of an expanding homoge- 
neous, isotropic general relativistic Friedman-Lemaitre cosmology. E. Lifshitz [3] found 
that disturbances with wavelengths larger than the 'dynamical' cosmological Jeans length 



grow like a power law, as compared to the background evolution, cf. [4,5]. 

The Jeans instability is generally regarded as forming the basis of our understand- 
ing of gravitational condensation. In particular, the non-relativistic criterion is invoked 
in astrophysical theories of the formation of stars and some smaller stellar systems, the 
relativistic one in cosmological theories about the formation of galaxies. 

Interestingly, while the analysis for the evolution of infinitesimal disturbances of a 
homogeneous and isotropic dynamical relativistic universe proceeds in an orderly manner 
(see in particular [4]), Jeans' original analysis for the evolution of infinitesimal disturbances 
of a homogeneous and isotropic static non-relativistic universe, which is reproduced in 
pertinent textbooks and monographs on the subject, e.g. [4-7], enjoys a rather questionable 
reputation. The following quotation is taken from the excellent monograph by James 
Binney and Scott Tremaine [7, p. 287ff.] (emphasis in the original): 

"We construct our fictitious infinite homogeneous equilibrium by perpetrating what we 
shall call the Jeans swindle after Sir James Jeans, who studied this problem in 1902 
(Jeans 1929). Mathematically, the difficulty we must overcome is that if the density and 
pressure of the medium Pq, Pq are constant, and the mean velocity Vq is zero, it follows 
from Euler's equation (5-8) that V$o ^ 0. On the other hand, Poisson's equation 
(5-9) requires that V $o ^ ^TlGpQ. These two requirements are inconsistent unless 
Po ^ 0. Physically, there are no pressure gradients in a homogeneous medium to 
balance gravitational attraction. A similar inconsistency arises in an infinite homoge- 
neous stellar system whose DP is independent of position. We remove the inconsistency 



by the ad hoc assumption that Poisson's equation describes only the relation between 
the perturbed density and the perturbed potential, while the unperturbed potential is 
zero. This assumption constitutes the Jeans swindle; it is a swindle, of course, because 
in general there is no formal justification for discarding the unperturbed gravitational 
field." 

Since modern cosmological applications do not any more rely on the original Jeans 
analysis but operate with the unproblematical relativistic dynamical analysis instead, and 
since a procedure for which "in general there is no formal justification" does not make 
much mathematical sense, one might want to conclude that the original Jeans analysis is 
best avoided at all. However, the importance of the stability analysis of an idealized infinite 
homogeneous and isotropic static universe lies not in real world physics applications but 
in the prospect of simplifying the analysis of certain physically relevant questions while 
retaining some essential features of a real system. As such it has been taught, tongue in 
cheek of course, to generations of students and newcomers to the field. But if we cannot 
backup the Jeans swindle by a formally correct analysis, then we face the dilemma of 
rendering a questionable service, surely to the novice in the field, but ultimately also to 
ourselves. 

As way out of this dilemma, Binney and Tremaine suggest that in special situations a 
formal justification of the "Jeans swindle" might exist. They list two examples [7, p. 288]: 

"However, there are circumstances in which the swindle is justified. For example, 
(i)... [if] ... the wavelength ... is much smaller than the scale over which the equilibrium 
density and pressure vary . . . the Jeans swindle should be valid for the analysis of small- 
scale instabilities. 

(ii) ... a uniformly rotating, homogeneous system ... can be in static equilibrium in 
the rotating frame and no Jeans swindle is necessary (although the stability properties 
are somewhat modified from those of the non-rotating medium because of Coriolis 
forces...)" 

Unfortunately, upon closer inspection neither point (i) nor (ii) really justifies the Jeans 
swindle. Thus, (i) would work if the effective Jeans length of some self-gravitating equi- 
librium were much smaller than the scale of non-uniformity of that equilibrium. However, 
the typical scale of non-uniformity of a self-gravitating equilibrium is precisely the effec- 
tive Jeans length, as emphasized in [2], so that any hypothetical "small scale instability" 
would have to be small in scale compared to the effective Jeans length, hence would not be 
analyzable by (1.1). In case of point (ii) the situation appears to be somewhat better, but 
also here there is a catch. In fact, while the introduction of uniform rotation with angular 
frequency vector Q does regularize the homogeneous gravitational problem in such a way 
that an analysis in the spirit of Jeans can be carried out without invoking any 'swindle' 
[8] , and while the resulting instability criterion for wave vectors satisfying k ■ S7 7^ is pre- 
cisely |kpCg — AnGpQ < 0, in agreement with (1.1), for this to justify the 'Jeans swindle' 
we would now have to be able to pass to the limit of a genuine 'Jeans-swindle-situation' 
from this 'no- Jeans-swindle-situation.' But exactly that is not possible, because the angu- 
lar frequency of a uniformly rotating equilibrium and the equilibrium mass density po are 
related by |Op = 2nGpo- As a consequence, the mass density po has to vanish in the non- 



rotating limit, and therefore a uniformly rotating system does not have a non-rotating limit 
in which the dispersion relation of the rotating system, which is different from (1.1) due 
to the presence of Coriolis forces [8] , would go over into the dispersion relation discovered 
by Jeans using his 'swindle.' 

In this paper we will present a mathematically clean vindication of the genuine "Jeans 
swindle," and a simple one at that, which emphasizes the physics underlying Jeans' rea- 
soning. We also mention a second, independent method which gives the same result. 

To begin with the physics, the first point to realize is that what counts dynamically 
are the forces, not the potentials. Hence, as long as we obtain a sensible dynamics in 
some sensible limit, we should not worry too much if some potential ceases to exist in 
the same limit. The second point to realize is that in Jeans' originally conceived system, 
three infinities are involved: (i) the infinite amount of matter, (ii) the infinite extend in 
space, and (iii) the infinite range of the Newtonian gravity. Of course, (i) and (ii) are 
not independent because the assumption that the putative equilibrium is homogeneous, 
and the perturbations only infinitesimal, couples these two infinities. However, (iii) is an 
independent source of trouble. Therefore, we are well advised to attempt the construction 
of the infinite system through a double limit. There are two mathematically natural 
possibilities: 

(A) We first let the mass and size go to infinity in a system with 'screened' gravitational 
forces, and subsequently 'switch on Newtonian gravity' by removing the screening. This 
treatment is borrowed from the class of infrared problems well known in quantum field 
theory. The standard procedure of handling infrared divergences is to apply an infrared 
regularization, to solve the regularized problem, and to remove the regularization at the 
end of the calculation, perhaps involving a 'renormalization.' As we will see in a moment, 
such a procedure works also here in an orderly manner. 

(B) We study the system from the beginning with classical gravity, though not in 
M? but in Sjj, which is §^ scaled to have radius R. Now the space is finite but without 
boundary, yet static. We have no problem in defining a homogeneous and isotropic static 
universe for classical gravity on S^j, and also not in studying the dynamics in its neighbor- 
hood. Letting R —^ oo subsequently, using an obvious renormalization, we arrive at the 
Jeans dispersion relation for a system in M.^, again in an orderly manner. 

In the rest of the paper, we will explain procedure (A) explicitly, leaving (B) as (easy) 
exercise for the interested reader. We first introduce the screened gravitational interactions 
and explain that this does not lead to physically unreasonable conclusions. Then we 
introduce the fiuid-dynamical equations of an asymptotically homogeneous system with 
gravitational interactions by taking the no-screening limit of the corresponding model 
with screened interactions. The linearized version is precisely the system discussed by 
Jeans, and his dispersion relation (1.1) follows at once. 

The equations of encounter-less stellar dynamics are treated analogously as limit of 
an asymptotically homogeneous system with screened gravitational interactions, yielding 
the corresponding Jeans dispersion relation for that model. 



2. SCREENED GRAVITATIONAL INTERACTIONS 

Consider first a sufficiently well behaved mass density function /9(x) which, for the 
moment, shall have finite mass, thus J^^ p{x) d x = M. Replace the familiar Newtonian 
potential of p at x G M"^, given by 

H^) = -G f ^-^p{y)d'y, (2.4) 

JRa |x-y| 

by the screened Newtonian potential 

v1>(x) = -g/ ^p{y)d'y (2.5) 

JRS |x-y| 

where k~^ is the screening length. In the limit of vanishing screening, i.e. k ^ 0, (2.5) 
reduces to the Newtonian potential (2.4). Therefore, if k is tiny enough, for instance 
K~^ = size of visible universe say, then for all concrete physical situations where one 
can apply the Newtonian gravitational potential (e.g., non-relativistic planetary motion, 
stellar dynamics and even some galactic dynamics) one can as well apply the above screened 
gravitational potential.^ 

Now drop the requirement that J pd x = M and consider a monotone sequence of 
mass densities that converges to a constant mass density po > (for instance, in sup norm, 
which means that sup^gRS \po — p(x)| -^ 0). Then (2.5) converges (likewise in sup norm) 
to a constant limit as well, given by 

*o = -Gpo / -. r- d^y = -AnGpo— (2.6) 

Jm^ |x-y| K^ 

while the Newtonian potential diverges, ^^ $ — *> — oo, as p ^ po- The divergence of $ as 
p ^ Po is not yet bad news, though, for we know what counts are not the potentials but 
the forces derived from them, viz. their gradients. Since a homogeneous system has no 
gradients, we would be in an acceptable limit situation if only we could guarantee that the 
gradient V$ would converge to zero as p — > po. Unfortunately, not only does the gradient 
V$ not converge to zero, its limit (whenever it exists) does not just depend on the limit 



^ Of course, this is not to say that a screened interaction with tiny screening enjoys an 
equal physical status as the Newtonian interaction - quite the contrary is true. By applying 
Occam's razor, there is no point in introducing screened interactions for the discussion of 
actual finite astrophysical gravitating systems in the non-relativistic regime. 

^^ In fact, $ = — oo even before p has converged to po, which is a consequence of 
considering convergence in sup norm and the fact that |x|~^ is not integrable over M^. If 
instead of sup norm convergence one uses the weaker concept of sup norm convergence on 
all compact sets (for instance by considering a mass density which equals po inside a ball of 
radius R centered at y and which vanishes outside the ball, and then letting R — > oo) then 
once again $ ^ — oo though this time $ stays finite for all finite R. Of course, \1/ ^ \l/o 
also now, this time on compact sets. 



density po but on the particular limiting sequence p ^ po (a fact which is well known 
in astrophysics and which has contributed to the general belief that Jeans' analysis is a 
'swindle.') 

The existence of the constant limiting screened potential \l/o on the other hand guaran- 
tees, in conjunction with the definition (2.5), that V\I/ -^ V\I/o = when p — > po- In other 
words, the screened gravitational interactions cancel themselves out when p = po. Since no 
pressure gradients are needed to counterbalance these self-balanced screened gravitational 
forces, it follows that such an infinite, homogeneous and isotropic fiuid is automatically in 
equilibrium. The infinite homogeneous and isotropic equilibrium fiuid with self-balanced 
gravitational forces can now be defined by simply taking the limit k ^ of this equilibrium 
family. 

In an analogous manner we can treat non-uniform mass density functions which differ 
from Po by the displacement of only a finite amount of mass. We could be more general, 
but this is certainly a reasonable class of mass densities to study. Hence, writing p(x) = 
Po + cr(x), the density disturbance o"(x) must be sufficiently integrable, satisfy 

a(x)d^x = 0, (2.7) 

and be bounded below by — po (for po + ^"(x) is a mass density and must therefore not 
be negative). For technical convenience, we actually demand that a be smooth and decay 
rapidly to zero at spatial infinity. 

The screened gravitational potential \1/ for such a mass density p(x) = po + o"(x) is 
readily computed. By the linearity of the integral formula (2.5), we have 

^(x) = ^o + V'(x) (2.8) 

where \l/o = —4'itGpo/k^ as before, and 

^{^) = -G f ^-^^a{y)d'y (2.9) 

Jr3 |x - y| 

Since a ^ constant, the gradient of \1/ in general now does not vanish, but is given by 

V*(x) = V^(x) (2.10) 

Unless the force density — p(x)V'i/'(x) is counterbalanced by pressure gradients, p(x) will 
not be a stationary density. The ensuing time-evolution will be studied in the next section 
using the conventional Euler, respectively Vlasov equations for such a system. 

The important point now is that because of the finite amount of mass involved in the 
density disturbance cr, the limit k ^ of Vifj exists and is given by 

lim V^(x) = V(/)(x) (2.11) 

where 

4>i^) = ~G [ ^—a{y)d'y (2.12) 

./in.3 x — y 



is precisely the expression for the Newtonian potential of the density disturbance invoked 
by Jeans. Therefore our nonlinear dynamical equations for a spatially asymptotically homo- 
geneous and isotropic fluid (respectively, stellar system) with screened interactions turn 
- without any swindle - into nonlinear dynamical equations for an infinitely extended, 
asymptotically homogeneous and isotropic fluid (stellar system) with Newtonian gravi- 
tational interactions. Their linearization gives the equations originally derived - in his 
heuristic manner - by Jeans. 

We conclude this section with a few remarks regarding Poisson's equation. We notice 
that the screened gravitational potential \1/ given in (2.5) satisfles the inhomogeneous 
Helmholtz equation 

A^ - K^^ = AnGp (2.13) 

Formally, as «: ^ 0, the Helmholtz equation (2.13) for \1/ goes over into the Poisson equation 

A$ = AnGp (2.14) 

for the gravitational potential $, provided the solution \1/ to Helmholtz's equation converges 
to a proper $ when k ^ 0. Now, for the asymptotically homogeneous system where 
P = Po + o- and \1/ = \l/o + ^, we have \l/o -^ — oo when k ^ 0, as obvious from (2.6). 
Hence, since \1/ y^ $ in this case, Helmholtz' equation (2.13) for \1/ does correspondingly 
not turn into Poisson's equation (2.14) for some $. However, once again this is not a 
problem because a constant potential is dynamically irrelevant. Since we are entitled to 
'renormalize' the potential by subtracting an overall (^-dependent) constant, we study 
\1/ — \l/o. The Helmholtz equation for \1/ — \l/o now does converge to a Poisson equation, 
namely the one which relates the density disturbance a"(x) to its gravitational potential 
disturbance (p. This is readily seen by recalling that \I/(x) — \l/o = i(^{^) is the screened 
gravitational potential disturbance associated with the density disturbance /9(x) — po = 
a"(x), so that the Helmholtz equation for \1/ — ^o simply reads 

Aijj - k'^i/j = ArrGa . (2.15) 

As K —* 0, t(j ^ (p , defined earlier in (2.9), and (2.15) turns into the Poisson equation 

A(j) = ATiGa (2.16) 

for the gravitational potential disturbance (f). All the ad hoc steps of the "Jeans swindle" 
have materialized in a mathematically clean way. 



3. THE FLUID-DYNAMICAL JEANS DISPERSION RELATION 

After our discussion of the screened gravitational interactions, the derivation of the 
Jeans dispersion relation (1.1) is now straightforward. The model considered in this section 
is the Euler model of the inviscid motion of a fluid with screened self-gravitation. The 
variables of the model are the fluid mass density p, pressure p, temperature T, fluid 
velocity u, and the screened gravitational potential \1/. The model equations comprise the 
continuity equation 

a,p + V-(pu) = 0, (3.17) 

the force balance equation 

p^tu + pu- Vu = -Vp-pV*, (3.18) 

an equation of state, for which (for simplicity) we choose the law of the classical perfect 
gas at constant temperature Tq, 

p = —pUbTq , (3.19) 

m 

and the Helmholtz equation for the screened gravitational potential \1/, 

A^ - K^-^ = AixGp (3.20) 

displayed here again for the sake of completeness. These equations have to be supplemented 
by asymptotic conditions at spatial inflnity, and initial conditions at some initial time, say 
to = 0. In particular, we demand that asymptotically at spatial inflnity the dynamical 
variables approach the values of the stationary, inflnite, homogeneous and isotropic equi- 
librium fluid in which the screened gravitational forces balance themselves. It is a trivial 
matter to verify that the set of constant variables, p(x) = po, p(x) = po — {po/m)kBTQ, 
u(x) = uo = 0, and \I/(x) = \l/o = —AttGpo/k^ for all x, forms such a stationary solution 
of (3.17-2.13). This constant solution will be our reference point. 

To inquire into the dynamics in the mathematical neighborhood of this constant solu- 
tion, we write p(x, t) = po + o"(x, t), and demand that the initial (t(x, 0) is smooth, rapidly 
decaying to zero at spatial inflnity, and satisfles 

(7(x,0)d^x = 0. (3.21) 

Then Jj^g (j(x, t)d x = for all t G (0, r), where r is the mathematical lifespan of the 
solution. Pressure and screened potential are represented accordingly, thus p(x, t) = po + 
a{'x,t)kBTQ/m, and\l/(x, t) = \l/o -|-'i/'(x, t). We also write u(x, t) = uo + v(x, t). Although 
for our choice of reference equilibrium we have uq = 0, whence u(x, t) = v(x, t), we prefer 
to introduce a new symbol for the deviation from the equilibrium velocity fleld simply as 
a reminder that more general equilibrium velocity flelds can be handled. 

Inserting the above representation of the dynamical variables into our fluid-dynamical 
equations, and already implementing our equation of state into the force balance equation, 
as well as using the fact that derivatives of constant functions vanish and that \l/o cancels 



versus po from Helmholtz's equation, we obtain the dynamical equations for the unknowns 

cr, V, V, 

dta + poV -v + V ■{av) = 0, (3.22) 

7 rri 

(po + a)dtv + {po + a)v-Vv = ^-^ Va - (po + ^) Vi/^ (3.23) 

m 

Atp - kV = 47rG'a . (3.24) 

All deviation variables are equipped with the asymptotic conditions that they vanish 
asymptotically as |x| ^ oo, allt. 

At this point already we can let k ^ in (3.22-3.24), thereby obtaining the nonlinear 
dynamical equations for the evolution of the disturbances of an infinitely extended fiuid 
with Newtonian gravity. The continuity equation remains unchanged, 

ata + poV-v + V-(av) = 0, (3.25) 

while force balance and Helmholtz equation become 

7 rj-i 

(po + a) atv + (po + a) V ■ V V = ^-^Va - (po + a)V(f) (3.26) 

m 

A(P = 4nGa. (3.27) 

All deviation variables are equipped with the asymptotic conditions that they vanish at 
|x| -^ oo. Notice that no linearization has been invoked so far. 
Proceeding on to the linearization of (3.25-3.27), we write 

a = (Ti+(T2 + ... (3.28) 

v = vi + V2 + ... (3.29) 

</> = (/>!+ </'2 + ... (3.30) 

where the index /c = 1, 2, 3, ... indicates the 'level of smallness.' Thus, a"2 is treated as one 
level smaller than ai; aiVcpi is at the same level of smallness as PQV(p2', etc. We are only 
interested in the first level of the hierarchy. Retaining only level 1 terms in (3.25-3.27), we 
obtain the dynamical equations at level 1, 

9tcTi + Po V • vi = , (3.31) 

7 rri 

Po 5tvi = -Vai - poV^i (3.32) 

m 

A(/)i = AnGai . (3.33) 

supplemented by initial conditions for ai , vi , and the asymptotic condition of vanishing at 
infinity for cri, vi, and (j)i. The solution of these equations is carried out in the standard 
way using Fourier and Laplace transforms, denoted by ^ and ~ , respectively. For the 
density perturbation in particular we find 

^ /, N c^ai(k,0)-pok.vi(k,0) 

ai{k,u) = -— ^ 3.34 

m ^/csTolkp — 47rGpo — a; 

from which we read off the original Jeans dispersion relation (1.1) for the disturbances of an 
infinite, homogeneous fiuid with isothermal equation of state and Newtonian gravitational 
interactions. We have done so without invoking a 'swindle,' or anything illegitimate of 
that sort. 

9 



4. THE STELLAR-DYNAMICAL JEANS DISPERSION RELATION 

In complete analogy, except that one has to be careful with the analytic continuations, 
one derives the stellar dynamical Jeans dispersion relation. It suffices to summarize the 
main steps. 

The dynamical variables of the model are, the particle density function /(x, v,t) on 
the one-particle phase space M"^ x M^ at time t G M, and the screened gravitational potential 
\&(x, t). They satisfy Vlasov's dynamical equations, which comprise the kinetic equation 

9t/ + V ■ V/ - V^ ■ d^f = (4.35) 

and the inhomogeneous Helmholtz equation for \1/, 

A^ - k'^^ ^ 4nG f f d^v (4.36) 



Our reference solution consists of a phase space density that is homogeneous in physical 
space with mass density po and Maxwellian in velocity space with temperature Tq, thus 
/ = /o(v) = /Oo(27r/cBTo/m)^/^ exp(— m|vp/2/cBTo), and the Helmholtz potential \1/ = 
\E'o = —AtvGpo/k'^, as before. Deviations from the stationary reference solution, written as 
/(x, V, t) = /o(v) + (/(x, V, t) and \I/(x, t) = \l/o + 'i/'(x, t), are required to approach reference 
values at spatial and velocital infinity. The dynamical equations for the unknowns g and 
t(j are 

O^g + ^.^g^^^. d^g = W ■ d^fo (4.37) 

A'iI;-k'^'i{j = 4tvG f gd^v (4.38) 

Taking the limit k ^ gives the nonlinear Vlasov-Poisson equations of an infinitely ex- 
tended, asymptotically (in space) uniform encounter-less stellar-dynamical system, 

dtg + v-Vg-V(P- d^g = Vcj) ■ d^fo (4.39) 

A(j)^47rG f gd^v (4.40) 

Expanding with respect to the levels of smallness, 

g = gi+g2 + ... (4.41) 

(I) = (pi + (p2 + - (4.42) 

and retaining only level 1 terms, we find the linearized Vlasov-Poisson equations, 

dtgi+y^-Vgi = V(Pi-d^fo (4.43) 

A(/)i = AtvG f gi d^v (4.44) 

10 



whose solution in terms of Fourier and Laplace transformation is standard. In particular, 
the phase space density perturbation, or rather its Fourier-Laplace transformed expression, 
for Xm{uj) < 0, reads 

«("• --') = , , ''^'V'?^ — rr^ r (4-«) 

which has to be analytically continued to Xm{uj) > 0. Apart from the ballistic term, absent 
in fluid theory, we immediately read off the stellar-dynamical Jeans dispersion relation for 
Xm{u) < 0, 

1 _ ±1^1 r -^^^ ^ (4.46) 

^ |k| Y 2fcsTo 

to be analytically continued to Im^u) > 0. In particular, using Plemilj's formula, we find 
that for Xm^uj) = the dispersion relation can be fulfilled only if 7^e(a;) = as well, which 
gives as stability boundary the familiar formula of Jeans, 



*^ (4,47) 

once again without any 'swindle.' 
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